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ABSTRACT 



Homogeneous Bethe-Salpeter equation for simplest Wick-Cutkosky model is studied in 
the case when the mass of the two-body system is more then the sum of constituent 
particles masses. It is shown that there is always a small attraction between the like- 
sigh charged particles as a pure relativistic effect. If the coupling constant exceeds some 
critical values there arise discrete levels. The situation here is analogous to the so-called 
"abnormal" solutions. The signature of the norm of these discrete states coincides with 
the "time-parity". The states with the negative norms can be excluded from the phys- 
ical sector-the one-time (quasipotential) wave-function corresponding to them vanishes 
identically. However the positive norm states survive and contribute to the total Green 
function (and the S-matrix) with the proper sign. 
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1 Introduction 



The main aim of this article is to demonstrate that equations of the relativistic quantum theory 
may have solutions which have no nonrelativstic analogue. Some of such solutions are well-known, 
e.g. the so-called abnormal solutions. Nakanishi wrote about the latters: "The non-relativistic 
common sense does not necessarily remain valid in the relativistic quantum field theory" 

The main ideas of this article were published as a preprint of Tbilisi Mathematical Institute in 
collaboration with A. Tavkhelidze and L. Vachnadze 0. 

At first the interest to the problem was initiated by the observed narrow peaks around 1.5MeV 
in the e+e~ distribution in the experiments of heavy ion collisions [Q . These resonances are known 
as the GSI resonances. They are situated above the 2me threshold and it naturally triggered 
interest to studies of bound states in continuum. In the papers of the Serpukhov and other groups 
1^ the single-photon exchange quasipotential equations in QED were studied and, relying on the 
numeric methods, large number of peaks were reported. 

The quasipotential equation can be obtained from the Bethe-Salpeter equation after certain 
approximation. So it was natural to try to cope with the same problem within the framework 
of the ladder approximation of the Bethe-Salpeter equation itself. It was particularly interesting 
because there was no such observation made on the earlier stages of studies of the Bethe-Salpeter 
(BS) equation (for a review see However, we think that nobody has searched for bound states 
in the continuum there. It is worth noting that when the total mass M > 2m the traditional 
method — Wick's rotation is, in general, inadmissible. So the problem must be considered in the 
Minkowski space. Structure of the ladder BS equation for this case was studied by Giinter and 
it was demonstrated that Wick's rotation is inadmissible. In the recent years interest to the BS 
equations in the Minkowski space has increased Below, employing standard methods, we 

will reduce the BS equation to the system of the one-dimensional equations and study unusual 
solutions in the continuous spectrum. Note that though further experiments have not confirmed 
existence of the GSI resonances the possibility of bound states embedded in the continuum is 
interesting on its own. 

2 Bethe-Salpeter equation in the Ladder approximation 
and its reduction 



We will consider the BS equation in the ladder approximation for a simple model 0,[llO| of two 
equal mass (m) scalar particles exchanging massless scalar meson, the so called Wick-Cutkosky 
model with the interaction Lagrangian: 

Lint = 5'l0102-4 + 5f20201^- 

Standard equation has the form (in the rest frame): 



V'(P) = ^ / d'^Q- 



P \^ 2 



2" — g j —m^ 



2 



Here P = pi+p2 is the total 4-momentum, p = ^{pi —P2) and q are relative 4-momenta, A = 
is the product of appropriate coupling constants and has dimension of the mass squared. In the 
rest frame where P = and pi = —p2 = p we will use the following notation for the total energy 
(mass) Pq = E = M = 2W. 
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Usually the Wick-Cutkosky model is considered as an auxiliary (fictitious) one. But in some 
reasonable approximation more realistic models reduce to the Wick-Cutkosky one. Indeed, in the 
single-photon (gluon) approximation in QED (QCD) for fermion and antifermion bound state BS 
amplitude we have an equation: 



where Xp = and the Feynman gauge is used. 

Usually a new function 



{p — g)^ + io 



1p{p) 



m + 7( p) 



is introduced M and the decomposition into the Dirac structures is considered: 



T 



Further one can obtain a chain of coupled equations for Dirac structures which can be con- 
siderably simplified by taking binding energies small (i.e. for weakly coupled system). All the 
components of \E' that survive in that limit of small binding energies (\l/^, \E'^) satisfy the 
same equations: 



[{p - qY + io] 



which coincides with the equations of Wick-Cutkosky model if we set AXprn? = A or Att = w?a 
where a = e'^/Air. Therefore for weekly coupled systems in single-photon (gluon) exchange ap- 
proximation in QED (QCD) we have in fact to deal with the equations of Wick-Cutkosky model. 

It is well known 0,0] that in the case W < m and A > (which corresponds to attraction) 
with small A eq. (|I|) has a discrete spectrum in accordance with the Balmer series for the hydrogen 
atom. Now our task is to show whether this equation has any discrete solutions when W > m i.e. 
M > 2m. 

It is convenient to transform the integral BS equation into a differential one in momentum 
space by using the identity: 



dpi dp^ 



1 



{p — qY + io 



We obtain: 



dpi 



-47i^iS^^\p-q). 



4Af*(p) 



(2) 



(3) 



{pI-P^ + W^- m2)2 - AplW^) 

This equation is mathematically equivalent to eq.(P if "^{p) satisfies the following boundary 
conditions |T^: p^\E'(p) remains finite for large values oi p and "^{p) must be finite when p = 0. 

We think that like in the Coulomb problem in our problem too it is possible to separate variables 
due to the hidden additional symmetry of the Coulombic interaction [jl3|. For earlier attempts of 



variable separation see [12 ,|13],[14 . 
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First of all let us single out the spherical angles in a standard way: 

vl/(p) =$(po,p.)l/™(^,0) 
where the spherical coordinates are defined as follows: 

Pi = Ps sin(6') cos(0), p2 = Ps sm{9) sin(0), p^ = ps cos{9), Ps 
Then eq.(^ is reduced to the partial differential equation with two variables: 



±\p\ 



Q2 Q2 

dpi dpi 



{PsHp) 



A\f{psHp) 



{pi — pI + — rn?y 
Boundary conditions for $(po,Ps) can be rewritten as 



AplW^) 



(4) 



Ps^{PO,Ps 



2 



Ps^{P0,Ps 



< oo 



PO^O 



p_ 

Ps 



iPs^{PO,Ps)) 



< OO 



(5) 



To investigate the energy spectrum in the regionl^ > m let us use the method of bipolar 
transformations , ll^ , [|14| . Peculiarity of our case is that the whole (po;Ps)-plane can be covered 



only by the means of several such transformations [|T4 

k sinh ai 



\^) Po cosh ai+cosh /3i ' P 

(2) V = fc sinh 02 „ 
V / Fv cosh a2— cosh /32 ' 



k sinh f3i 



cosh a 1 +cosh f3i 
k sinh 132 



|a2| 7^ \P2\ 



( 'Q\ ^ _ k cosh 03 
^ ^ F^ sinh rvQ-t-sin}' 

(4) Po 



sinho3+sinh/33' 
k cosh a4 



Ps 
Ps 



cosh 02 —cosh (32 

as ^ -(3; 



sinh a3+sinh f3'i 
k cosh /34 



^4 



3 

/?4. 



sinh a4+sinh /34 ' sinha4+sinh/34 

The hyperbolic angles in all these transformations run over the whole infinite intervals: 

— oo < ctj, Pi < oo 

while the parameter k is determined by the condition of separabilty in eq.(0) and equals to {W"^ — 

Action of these transformations is shown in Fig.l. After performing these transformations 
variables in eq.(^) separate: 

Ps^{Pq,Ps) = f{ad9{(^i) « = 1,2,3,4 
and we obtain the following equations for each of the functions / and g: 
(a) For transformations (1) — (2) (the shaded regions in Fig.l): 



da-' 



+ hi- 



sinh^ /3 



a — cosh a 



9{(3) = 0, 
/(«) = 



0; 



(6) 
(7) 
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(b) For transformations (3) — (4) (the unshaded regions in Fig.l): 

g{P) = 0, 



+ /i2 + 



+ h2 



cosh^ P 

\/m? 
a + sinh^ a 



/(«) = 0. 



(9) 



In these equations 



> 1 



and hi , h2 are the separation parameters (new quantum numbers) . The fact of separabihty is 
provided by the symmetry of relativistic problem in the 4-dimensional world |T^. 

Equations obtained above describe one- dimensional motions in the space of hyperbolic angles. 
Equations for g{j3) are kinematic ones defining the values of separation parameters hi 2. Equations 
for f{a) contain dynamical quantities (A, a) and hence must allow to find the energy spectrum. 

Boundary conditions for g{(3) and f{a) follow immediately from (^ and are the same in different 
regions of the {po,ps)-pla.iae 



(7(0) = (7(±oo)<oo 
/(0)<oo /(±oo) = 



(10) 
(11) 



Note that the boundaries of different regions in (po,Ps)-plane can be approached by taking 
limits a,f3 ±00 in the relevant transformations (1) — (4). 



3 Solutions of the one-dimensional equations 

First consider the pair (|D-(0). Eq.(||) with the boundary conditions ( |10|) has the following solution: 

g{f3) = gi{f3) ^ Lmh^^' cosHVhP), h>0 

I [12) 
0, hi<0 

In general, the second equation of this pair can give the relation between energy W and the 
parameter hi 

W = pnihi), n = 0,1,2... 

Despite the fact that n can take discrete values, energy is not quantized because hi remains to 
be an artbitrary positive number. 

Things are essentially different in the unshaded regions (3) — (4). Eq.(^ has solutions with 
positive and negative signs of /i2- The positive h2 are not resulted while the negative ones are 
quanttizes for / > 0. 

h2 = -{l-Nf, N<1. 

where N takes odd integer values. 

Next consider eq.(^ in more detail. Expression 



a + sinh a 
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plays the role of "potential". 

For positive A (i.e. attraction) V{a) is the positive definite and convex function of a (Fig. 2). 
Obviously we can have only continuous spectrum in this case witth posititve h2- 

When we take A negative (repulsion) V{a) becomes negative definitte and concave function of 
a (Fig. 2). Therefore in this case there appears principal possibility to have discrete eigenvalues 
provided the depth of V{a) is high enough for a given negative value of h2 = —{I — . 

One can examine easily that continuity on the boundaries of shaded and unshaded regions 
imposes no new constraints on separation parameters hi 2- So there is no discrete specttrum for 
positive values of A when W > m. 

Now we turn back to eqs.(H)-(|^) with /i2 < and A < 0. Eq.(|D is a Heun's differential equation 



15]. Its solution with h2 = {I — N)"^ and A = — |A| < imposed by boundary and normalization 
conditions leads to complicated transcendential equations. To estimate qualitatively the coupling 
constant dependence of the energy eigenvalues one can make use of tthe following approximate 



expression 



As > it follows that 



1) 



z/=l,2, 



\2 

> z/(z/+l) =2,6,12, 



m? 



i.e. discrete levels apeear when |A| exceeds some critical values. It is worth noting that for some 
fixed values of A^/m^ there is finite number of discrete levels. For example, when 



there is only a single level. And for 



2<^<6 



A2 

6 < — < 12 



the number of levels equals 2, etc. 

If we take |A| in accordance with Balmer formula mentioned above in the following form Avr = 
m^Z^a, then the above evaluation can be cast into 

> u{u + 1) 

TT 

i.e. Zrriin > 30. 

We want to note tthat eqs.(||)-(P) were derived earlier in [|l^] but with erroneous sign of A. 
Hence author deduced bound states for positive values of A. 



4 Normalization condition 

Let us now discuss the normalization condition for the BS amplitude, which is nonvanishing only 
in the unshaded regions and has the form: 

,/ N {smha + smhpYf{a)g{f3) 

"f^P) ~ _ ^2)3/2 cosh=^ ^(^2 + ^2 gi^h^ «) ^'^^ 
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The role of normalization is to provide the correct relation between the aplitude ifj and four- 
point Green function. The exact form of the normalization condititon for the model under con- 
sideration looks like 

z I d'pmGt{p)m = ^b\^^^ (14) 

where ^^(A) = iW^^X) and is the so called norm factor determining the sign of the pole 
contribution in the total Green function at S* ~ Sb{X) 

G ~ ieb — 5 5 + Reg (15) 

O — OB 

Integration in the normalization condition (|l^ is carried over the unshaded region only. 
The spectral representation for ip{p) and ip{p) that we are going to use has the form |^: 

'^{P0,P) = -TT- / dqo — + — / dqo -; 16 

2m Jw+{p) Po — Qo + 2m J-oo Po — Qo — io 

^{Pq,P) = -TT- / dqo — + — / dqo 17 

2m Jw+{p) Po — qo + io 2m J-oo Po — % — io 

where w^{p) = W — {p^ +m^)^/^, w+{p) = —w^{p) are the enndpoints where cuts in the complex 
Po-plane extend from. As long as 3-momentum p satisfies < < — these two cuts overlap 
(Fig.3). 

The spectral reprezentations (p!6D-(p!7D imply 

ij{z) = -r{z*), z = npo + i^po (18) 

There are "dangerous" regions in the normalization integral where left-hand and right-hand 
cuts overlap. This overlapping happens when 

Po e —^0,0 for right-hand cut, 
Po & 0,wo,0 for left-hand cut. 

In these intervals < p^ < k"^ = W"^ — m? . In the (po)Ps)-plane "dangerous" points are located 
between the profiles of functions w+{p) and W-{p) for pi < k"^. This region is drawn in Fig.4. 

Note that the "dangerous" points completely lie in the region where our BS amplitude vanishes 
identically. Therefore, owing to this specific behaviour of the BS amplitude for bound states with 
W > m, intervals of po where cuts overlap do not contribute to the normalization integral. 

It seems to us that this property is very important and can show up in more interesting cases 
too. Specifically this fact may be employed for confining kernels, because sometimes bound states 
above threshold {W > m) do arise for such kernels. Example considered above indicates that 
in similar cases the (po;Ps)-plane may bear some structure — namely the BS amplitude may 
have support region in the (po,Ps)-plane. So effectively we can take that cuts in the spectral 
representation start from the origin {po = 0) and, consequently, perform Wick rotation without 
any complications. This is the case when Wick rotation is admissible for dynamical reasons. As a 
result we obtain: 

d^pdpi%l){ipi,p)GQ^{ipi,p)i){ip^,p) = e\—^^ (19) 
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Now, according to (|T8|) and symmetry of eq.(^ under the relative time reflection po ~^ ~Po we 
have: 

i^{m,p) = i^*{-ip4,p) = -nBV^*(ip4,p), (20) 

where denotes the so-called relative-time parity ( "po-parity" ) of the bound state. 
Using (pOD in the normalization integral we get 

d^pdp4'ijj{ip4,p)GQ'^{ip4,p)'ip{ip4,P) = 

As long as GQ^iip^.p) is positive definite in the unshaded regions the full integral above is positive 
definite too. Besides, from the explicit form of our solution we know that 

and therefore with appropriate choice of overall multiplicative constant up to which the BS am- 
plitude is determined the norm factor sb coinsides with "po-parity" of the bound state 



B 



In other words the states with positive (negative) "po-parity" have the positive (negative) norm. 

States with the negative norm can be be eliminated in from the physical sector as the corre- 
spondinng "one-time" quasipotential wave functions vanish identically: 

dpoipiPo,P} = . 

On the other hand the positive norm states, i.e. the states with positive "po-parity" do survive 
and produce the pole contributions to the total Green function (and S-matrix) with the correct 
sign. 



5 Conclusions 

We have demonstrated that in the BS equation for the equal mass scalar particles' repulsive 
(A < 0) interaction via scalar massless particle exchange there emerges an effective attraction 
above thershold. This attraction can lead to existence of bound states for values of the coupling 
A larger then the critical one: 

|Ac| 



> 2 . 

In QED with the one-photon exchange this corresponds to 



^>2. 

TT 



It is worth mentioning that attraction is present for all values of A as a purely relativistic effect 
in the kinematic region where the retardation effects can not be treated by perturbative methods. 
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Figure 1: Regions in Pq, Ps plane corresponding to bipolar trans- 
formations (1-4) 
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Figure 2: The "potential" V{a) for both signes of A. 
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Figure 3: Complex pq plane. Here ujq = [u!-{p)max = W — m > 0. 
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Figure 4: "Dangerous" region in the po, Ps plane. 
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